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Abstract 

Discrete-time quantum walk evolve by a unitary operator which involves two operators 
a conditional shift in position space and a coin operator. This operator entangles the coin 
and position degrees of freedom of the walker. In this paper, we investigate the asymptotic 
behavior of the coin position entanglement (CPE) for an inhomogeneous quantum walk 
which determined by two orthogonal matrices in one-dimensional lattice. Free parameters 
of coin operator together provide many conditions under which a measurement perform 
on the coin state yield the value of entanglement on the resulting position quantum state. 
We study the problem analytically for all values that two free parameters of coin operator 
can take and the conditions under which entanglement becomes maximal are sought. 
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1 Introduction 

Among classical algorithms, many are based on classical random walk. Markov chain simula- 
tion, which has emerged as a powerful algorithmic tool [1] is one such example. Like classical 
random walk, the quantum version of it has also become an important constituent of quantum 
algorithms and computation. Though quantum random walk was first introduced by Aharonov 
et al. [2], but unlike its classical counterpart, the evolution of the quantum version is unitary, 
reversible and has no randomness associated with it during the evolution. Therefore, keeping 
away the term 'random', quantum walk has been the preferred usage. In the classical random 
walk the particle moves in the position space with a certain probability, whereas the quantum 
walk, which involves a superposition of states, moves by exploring multiple possible paths si- 
multaneously with the amplitudes corresponding to different paths interfering. This properties 
makes the variance of the quantum walk on a line to grow quadratically with the number of 
steps (time), compared to the linear growth for the classical random walk [3, 4]. Quantum walk 
is studied in two standard forms: continuous-time quantum walk and discrete-time quantum 
walk that was introduced by Farhi and Gutmann [5] and Watrous [6], respectively. In the 
continuous-time quantum walk, one can directly define the walk on the position Hilbert space 
[5], but in the discrete-time quantum walk, in addition to it is necessary to introduce a coin 
Hilbert space to define the direction in which the particle has to move [3]. Due to the coin de- 
gree of freedom, the discrete-time variant is shown to be more powerful than the other in some 
contexts [7]. Childs describe a precise correspondence between continuous- and discrete-time 
quantum walks on arbitrary graphs [8]. Both these types of quantum walk have been widely 
used in algorithms for a diverse of problems. See for example [9, 10, 11, 12, 13, 14]. Beyond ap- 
plications in quantum algorithms, quantum walk is emerging as a potential tool to understand 
various phenomenon in physical systems and has been employed to demonstrate coherent con- 
trol over quantum many body systems. See for example [15, 16, 17, 18]. Some experimental 
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progress on the implementation of quantum walk has been reported [19, 20, 21, 22, 23, 24]. 
Entanglement is one of the attractive properties of quantum that does not appear in classic 
and is a very useful resource to perform various quantum tasks, a recent review by Horodecki 
et al. [25] discuss many of these aspects. Unfortunately, physical limitations and noise effects 
change the amount of entanglement and restrict it's efficiency One way of circumventing this 
problem would be to generate entanglement. Quantum walk is one such process in which 
the conditional shift operator, which governs the itinerary of the quantum walker, induces 
entanglement between the degree of freedom of the coin and the spatial degree of freedom of 
the walker. This entanglement fluctuates with each step and eventually settles down to an 
asymptotic value that depends on the initial state of the quantum walk and on any bias in 
the quantum coin operator. Entanglement is a basic resource in quantum algorithms, and 
further work is required in order to fully understand its properties in the QW. While, entan- 
glement between coin and position of the walker, as a type of quantum correlation, distort 
the distribution of the quantum walk and give it peaks and troughs, especially at the ends 
of the top hat [31, 32], it is necessary to obtain linear spreading and mixing times [33] and 
can be use to gauge the impact of the added randomness and decoherence [30]. In addition, 
it can be exploited for quantum information and communication purposes [26, 27, 28, 29]. 
Entanglement has also given us new insights for understanding many physical phenomena in- 
cluding super-radiance [34], superconductivity [35], disordered systems [36] and emerging of 
classicality [37]. In particular, understanding the role of entanglement in the existing methods 
of simulations of quantum spin systems allowed for significant improvement of the methods, 
as well as understanding their limitations [38]. So, studying entanglement during the quantum 
walk process will be useful from a quantum information theory perspective and also contribute 
to further investigation of the practical applications of the quantum walk. Carneiro et al. 
[41] studied the long-time asymptotic coin-position entanglemant of quantum walks on various 
graphs (Z , Z 2 , triangular lattices, cycles). Venegas-Andraca and Bose [42] also investigated 
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generation of entanglement between two walkers. Using Fourier analysis techniques, Abal et 
al. [43] analytically computed asymptotic coin position entanglement of the Hadamard walk 
on one dimension for both localized and non-localized initial conditions and in the same way 
Annabestani et al. [44] studied asymptotic entanglement in a two-dimensional quantum walk. 
Chandrashekar et al. [45] consider a multipartite quantum walk on a one-dimensional lat- 
tice and studeid the evolution of spatial entanglement, entanglement between different lattice 
points. 

Most of the original work considered homogeneous walks, which the amplitude for moving does 
not depend on position. The idea of looking at inhomogeneous walks is not new, the recognition 
that it is natural to allow coins to be position dependent may be found in [46, 47, 48, 49, 50, 55]. 
To probe the possible long time behaviours and quantify the entanglement between the coin 
and position of such walks is our motivation for this work. In this paper, we investigate the 
asymptotic behavior of the CPE for an inhomogeneous QW which is determined by two orthog- 
onal matrices in one-dimensional lattice. Our inhomogeneous QW is a kind of the generalized 
model defined in [55, 56] which is based on the idea of the Aubry- Andre model [57]. The limit 
distribution of the this QW has the probability density given by Dirac delta function which 
is called a localization for the inhomogeneous QW (for more details see Ref.[55]). We show 
how the entanglement approaches asymptotic values depending on the choice of initial state 
and coin bias with study the problem analytically for all values of 9 and 9\ for large number 
of QW steps. We find that CPE is different for odd or even positions and when on the coin 
parameters took the precise value | walk alter to a bounded motion. 

The plan of the paper is as follows: in section 2 the model for discrete-time quantum walk is 
presented in detail. Section 3 defines the entanglement and provides the required formalism 
in the Fourier space leading to the long-time reduced density operator for arbitrary coin op- 
erations and initial states. Our inhomogeneous walk with two-period has defined in section 4 
and the asymptotic CPE entanglement of our walk as a function of coin operators for local 
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and non-local initial conditions is obtained. Conclusions are given in the last part, section 5. 

2 Definition of the quantum walk 

The discrete-time quantum walk (DTQW) on a line is defined on a Hilbert space H = H p cg>H c , 
where H c is the coin Hilbert space, spanned by the basis state \L) = [1,0] T and \R) = [0, 1] T , 
where T denotes the transposed operator, that represents two sides of the directions of the 
motion and H p is the position Hilbert space, spanned by the infinite basis states {\x); x e Z} 
that represent the position of the walker. In one step of the QW we first make superposition 
on the coin space with a coin operator C 6 17(2) and after that we move the particle according 
to the coin state with the translation operator S as U — S.(I P <g> C), where S is defined as 
S = T, x \x — l)(x\ <g) \L)(L\ + \x + l)(x\ ® \R){R\, and I p is the identity operator in H p . The 
most widely studied form of the DTQW is the walk using the Hadamard operation as quantum 
coin. It is an unbiased coin operation, the resulting walk is known as the Hadamard walk, 
which determined by the Hadamard gate: H — ^ 
The state of QW at time t is expressed by 

!*(*)> = EI*>®WM)>, (2- 1 ) 

X 

where \ip(x,t)) = a(x,t)\L) + b(x,t)\R) denotes the coin state and a(x,t) (b(x,t)) is the am- 
plitude of the base \x, L) (\x, R)) at time t which belong to the complex number C, satisfying 
the normalization condition J2 X \ a { x ^)\ 2 + \b(x,t)\ 2 = 1. Total state after one steps, is given 
by \^(t + 1)) = U\^(t)). The Fourier transform, as first noted in this context by Nayak and 
Vishwanath [4], is extremely useful when single-step displacements are involved because the 
evolution operator is diagonal in /c-space. The Fourier transform \ip(k,t)) (k e (— 7r, 7r)) is 
given by 

|^(M))=5>-^(*,t)>, (2-2) 

X 
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and by the inverse Fourier transform, we have 



^-e ik ^(k,t)) 



(2-3) 



The time evolution of \ip(k,t)) after one step of the walk is \tp(k,t + 1)) = U\ip(k,t)) , where 




, e ik o 

U = R(k)U and R(k) = 

V e" ! 



3 Entropy of entanglement 



As mentioned above, the conditional shift in evolution operator of a DTQW entangles the 
coin and position degrees of freedom of the walker. To efficiently make use of entanglement 
as a physical resource, the amount of entanglement in a given system has to be quantified. 
Entanglement in a pure bipartite system is quantified using standard measures known as 
entropy of entanglement corresponds to the von Neumann entropy [51]. It is a functional of 
the eigenvalues of the reduced density matrix and is given by the formula: 



In this equation, p c = tr p (p) is the reduced density operator obtained from p = (UYp(t = 
0)(Wy by tracing out the position degrees of freedom. Since p c has two dimension, this 
quantity is Se £ [0, 1], i.e., Se = for a product state and Se — 1 for a maximally entangled 
state. Note that, in general tr(p c ) = 1 and tr(p 2 c ) < 1. The entropy of entanglement can 
be obtained after digitalization of p c . This operator which acts in H c is represented by the 
Hermitian matrix as 



tr(p c log 2 p c ). 



(3-4) 




(3-5) 



where 



a(t) =^|a(x,t)| 2 



^ dk 



a(k,t)\ 2 , 



X 
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7 (t) = X>(M)l 2 = £^(M)I 2 - 

Where * denotes the conjugate value. The eigenvalues ri )2 of p c , can be computed as 

r lf2 = |[1 ± ^1 + 4(|/3(t)P - «(t) T (t))]. (3-6) 

Therefore, by using Eq. (3-4) one can obtain the entropy of entanglement as 

S E = -(rilogri + r 2 logr 2 ). (3-7) 

To study asymptotic entanglement for Hadamard walk, this method has been also used by 
Abal et al. [43]. They show analytically that for localized particle with every initial coin 
states, the asymptotic entanglement is 0.872 and when nonlocal initial conditions |\l/±(0)) = 
|(| — 1) ± |1)) eg) (\L) + i\R)) are considered, the asymptotic entanglement varies smoothly 
between almost complete entanglement, Sg 0.979, and no entanglement (product state), 
« 0.661. 

4 Two-period QWs 

When one quantum coin acts for all the walker position, we face with a homogeneous walk 
which the amplitude for moving does not depend on position but inhomogeneous walk differ 
from it, in that we allow the coin operator to depend on the walker position. Now we will 
define an inhomogeneous quantum walk in a similar manner to the standard quantum walk, 
which the coin operator C to be dependent on x. In the case of walks on the line, C x could 
be an arbitrary unitary operator on the two-dimensional coin space. Indeed, there is no need 
to restrict to walks in which only moves to the nearest neighbours, more generally one could 
allow there to be transitions from a given point to any other point on the line. However for 
the purposes of most of our discussion of walks on the line we will focus on the simplest case 
of a transitions to nearest neighbors. Let {C x ;x e Z} be a sequence of orthogonal matrices 
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with C2S = Hq and C2.S+1 = Hi (s e Z), where 

(cos 6* 7 sin # 7 
I , (4-8) 
sin # 7 — cos # 7 

which 7 = 0, 1; O — 27r£(2s), ^ = 27r£(2s + 1) are free parameters of coin operators which can 
be altered to choose the quantum coin operation and ( G R is the inverse period of the coin 
operations. With this selection of coin operators, we face with a two-period quantum walks 
which are determined by two orthogonal matrices, and are also self-dual under the Aubry- 
Andre duality. The inhomogeneous quantum walk is restricted to the finite interval [— Q,Q], 
when C = 47? with relatively prime P (odd integer) and Q. So, the limit distribution of the 
this inhomogeneous quantum walk, divided by any power of the time variable is localized at 
the origin [55]. Limit probability distribution of this two-period QWs computed by Machida 
et al. [50]. In the rest of this work, we shall be concerned with clarifying the asymptotic value 
of S e for both local and nonlocal initial conditions. 

4.1 Asymptotic entanglement from local initial conditions 

Let us first consider in detail the simple case of an initial state localized at the origin with no 
CPE as |^(0)) = |0)<8#(0,0)), where |^(0,0)) = a(0,0)|L)+6(0,0)|i2) is the coin state. Below, 
we provide an analytical explanation for the observed values of asymptotic entanglement in 
the local case. The time evolution for this initial state in /c-space regarded as 

\i>(k, 2t)) = (H^oYltpik, 0)) for even times 

(4-9) 

\i/>(k, 2t + 1)) = HofaHoYMk, 0)) for odd times. 
We note that 0)) = \tp(0, 0)) and H 1 = R(k)Hj. Two eigenvalues of HxHq are given by 



A 7 (/c) = cqC\ cos2A; + s si + (— lY'iy 1 — (c ci cos2A; + so s i) 2 for (7 = 0, 1), (4-10) 

where c 7 = cos# 7 and s 7 = sin# 7 . The eigenvectors |V^(fc)) corresponding to X 1 (k) are 

/ u{k) \ 

|v;( ^H W + m)mJ' (4 ' u) 
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which the elements of this matrix are as follows 

u(k) = SQC\e 2lk — cqS\ 

v{k) = —icqCi sin2/c 

w(k) = i\Jl — (c Ci cos 2k + s Si) 2 , 

and iV 7 is the normalization constant. Using the spectral decomposition for H\Hq, the Fourier 
transform \ip(k,2t)) is expressed as 

\4>(k,2t)) = j2\i f (k)(V,(k)$(k,0))\V,(k)), (4-12) 

7=0 



and its spinor components are obtain 



a(k,2t)= u mF(k) + ^G(k)), 



(4-13) 



b(k, 2t) = v(k)(^F(k) + X -^G(k)) + w(k)(^F(k) - ^G(k)), 
where 

F(k) = u*(k)a(k, 0) + (v*(k) + w*{k))b{k, 0), 
G(k) = u*{k)a{k, 0) + (v*(k) - w*(k))~b(k, 0). 
Quantum walk after 2t + 1 steps is 

|^(fc,2t + l)> = 2^(fc)(V 7 (fc)|^(fc,0)>ffo|V 7 (*)>. (4-14) 

7=0 

It is easy to verify that coefficients of \L) and \R) have the forms as 

a(k, 2t + 1) = e ik {(c u(k) + s v(k))(^F(k) + ^G(k)) + s w(k)(^F(k) - ^G(k))} 
b(k, 2t + 1) = e-^{(s u(k) - c v(k))(^F(k) + ^G(k)) - c w(k)(^F(k) - ^G(k))}. 

The relevant quantities for the entropy entanglement are a(t) and (3(t) which defined in Eq. 
(3-5). , after 2t steps we obtain 
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«(2t)=^f{KA:)| 2 (^f + 

+Sf(A 2 1 t (fc)F(A;)G'*(fc) + \?(k)F*(k)G(k))}, 

(J(2t) = %{u{k)v*{h){\^f + ) + «(A>-(*)(1^ - Mf) 

+ ! fSr 1 ( A o t (^)^(A:)G*(A;) + Af(fc)^(fc)G(fc)) - ^^(Xf(k)F(k)G*(k) - Xf(k)F*(k)G(k))}. 

(4-16) 

The time dependence of these expressions vanishes in the long time limit by using Riemann- 
Lebesgue lemma (see Appendix), and the asymptotic values for even time steps are 



« = /-.fKA:)| 2 (^f + ^), 
? = §Mk)v*(k)0-Zf + ) + «(JbX(Jb)(J^l! - )}. 



(4-17) 



Where we use overline to indicate that the asymptotic limit has been taken, i.e. ~p c = lim^oo p c . 
After some algebra and use the Riemann-Lebesgue lemma the asymptotic values a and (3 for 
odd time steps can be obtained as 

« = g {(c 2 \u(k)\ 2 + c s (u(k)v*(k) + u*(k)v(k)) + st(Wk)f + HA;)| 2 ))(^ + J^) 

+(coS (u(k)w*(k) + u*(k)w(k)) - 24v(k)w{k))(Mf - J^)}, 
P = S% ^ k {(c s (\u(k)\ 2 - \v(k)\* - \w(k)\ 2 ) - clu(k)v*(k) + slu*(k)v(k))(^ + ) 
+(-c§u(A>*(*0 + s 2 u*(k)w(k) + 2 W (*>(*0)( J£ P ! - 

(4-18) 

In order to further illustrate the effects of bias coin on the asymptotic entanglement level, the 
asymptotic entanglement corresponding to initial condition given by |^(0)) = |0) ® ^^/^ is 
shown in Fig.l for even times and Fig. 2 for odd times, discretely. The plots in these figures 
show the entropy of entanglement as a function of 9q and 9\ defined in Eq. (3-7). As it is 
tangible from the form of coin operator, the period of asymptotic entanglement variations is -n 
and the mirroring about the vertical axis indicate that Se(0 o , 6\) is an even function. Fig.l (A 
and B) shows that, at the exact value 9 1 = n/2 and all possible values of 9 since a — \ and 
(3 = ^ so the eigenvalues of p c are r± — 1 and r 2 = 0, thus we face with a completely separable 
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Figure 1: (Colour online) Asymptotic of entanglement Se as a function of 6 and 9\ for 
ao(k) = ibo(k) = ^= after even time steps. 

state and entanglement drops to zero, S^ ter even ttme steps \6 , f ) = 0. Moreover for 6 = n/2 
and full of Q\ we obtain a = \ and ft = -r thereupon r± — | and r 2 = ~, in this case we 
have 4 a/ter even time steps \l,6 l ) = 0.811278. In comparison, we find that, S% fter odd Ume steps) 
versus 6q an d $i is symmetric and when these variables take the accurate quantity |, we have 
a. = \ and (3 = therefore maximal entanglement shows itself, see Fig. 2 (A and B). To 

elaborate this further, we will consider the steps of evolution of the DTQW using quantum 

/0 1\ /cos#! sin#! \ 

coin operations given by H = and H\ = . After several steps of 

\1 0/ \sin0i -cosOj 

the DTQW, the state can be written as 

[V(2*)> = ^{(e^ 1 + (-l)V iei )|0) (8) (\L)+i\R)) 
+( e -«ffi - (_i)* e «tfi)(| _ 2) <g> |L) + z|2) ® 

(4-19) 

| V(2* + 1)> = ^{(e^ 1 + (-l)V tei )(^| - 1) ® \L) + |1) <8) |it:)) 
+ ( e -«* _ (_i)V* ei )(«|l) <8) |L) + I - 1) ® \R))}. 

We can observe that, when H reduces to the Pauli X operator, then the inverse period of the 
coin operations corresponds to so our inhomogeneous quantum random walk is restricted 
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Figure 2: (Colour online) Asymptotic of entanglement Se as a function of #0 and Q\ for 
aa(k) = ibo(k) = after odd time steps. 

to the finite interval [—2,2] [55]. The reduced density operator can be computed from above 
equation: 



p c (2t) 



;i + (-i^cos^o) - 



+ cos{2t9 1 )) 



and p c (2t + 1) 



\ 
j 



after averaging over t in the long time, the matrix element of the reduced density operator is 
obtained: 



p c (after enen steps) 



2 4 



i I 
4 2 



and p c (after odd steps) 



\ 



It is easy to check that S^ teT even Ume steps) (f , 6 l ) = 0.811278 and S^ fter odd tme steps \^ 6,) = 
1, see Figs.l and 2. When we let Hi diminish to coin shift operator and other parameter of 
coin (#0) alter, the wave functions of the QW can be written as 

M2t)) = ^J*>\0)®(\L)+i\R)) 

Vi (4-20) 

1-0(2* + 1)) = i=*£ e <(*+i)«°(| - l) ® \L) - i\l) <g> |it!)). 
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It is discernible that particle remain in the origin in all even times as a separable state and in 
all odd times, it with maximal entanglement exist in the positions —1 and +1 with coin states 
\L) and \R), respectively. Hence, the quantum walk is bounded. 

One of the most remarkable facets of this walk, is behavior of asymptotic CPE for — O < 
6i < O in all even times. In this interval S^ ter even Ume steps ^ gets stuck value which is 
determine by biases of Hq, this characteristic ascertain that when 9q attains constant value, 
conversions of X can not change the level of entanglement. The change of asymptotic CPE 
is from X G (0 , 7r — O ), in this interval S^ ter even Ume steps ^ fluctuate between the maximum 
and minimum possible entanglement levels and its quantity depend on both parameters of 
coin operator, O and 0i (see Fig. 1 (A)). Moreover, with increase the value of O from zero 
to |, maximum quantity of long time limit entanglement reduce. For example we study the 
behavior of asymptotic entanglement after even time steps for O = j and we find that a = \ 
and for 9\ = ~^f, ^p, f3 = |(1 — \/2) but against other values of 9\ exclusion of ^ we have 

7j = J- S ec 3 6>i {3«7T cos 9 X + in cos 30 x + 2 cos 5/2 



3 cos 9 1 - cos 39, + VI + cos 29, sin 6,{ln\ , m^^^^^ M \ 

\/ 6 cos 2 0i — 2 cos 0\ cos 30i +4 sin 26Wl+cos 26»i ^ 
_i r -t(l+cos20i-sin20i-2Vl+cos20i) i-| 

y/e cos 2 0i - 2 cos 6»i cos 30i +4 sin 20 x ^1+cos 20i 



+1 J -3 COS 9 1 + COS 3fli + Vl + COS 20! Sin 9~{\n[ . -^(l+c°s20l-sin20i+2Vl+cos20i^ 
v V - 3 cos 2 6»i +cos 6»i cos 30i +2 sin 26»i V 1 +cos 26»i 

_ ^ ^(1+008 201-8111201+2^1+008 201) jj^j 

\f— 3 cos 2 0i +cos 0i cos 30i+2 sin 20i -v/l+cos 20i ^ 



when < 0x < f , the exact value of asymptotic CPE is S% jter even Ume steps \%9 l < f) = 
0.872429, this precise value are coincident with that obtained for homogeneous Hadamard 
walk with local initial condition [43], from | to ir — | asymptotic CPE variate slowly and 
achieve zero for correct 9, = While 9\ close to n — |, it gradually augment to 0.872429 
when the parameter 9\ attains n — | and the asymptotic entanglement maintain this value to 
0i = n. In plot (B) presented in Fig. 1, overlapping graphs endorse our notion that for 0i < O 
the limiting value of the entanglement is regulated with coin operator H . For example after 
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even time steps we obtain ol—\ and for 9 = ^p, ^p, (3 = |(1 — V2) but versus other values 
of #0 exclusion of ^ we have 

P = — 1 ={i sec 3 / 2 6W-3 cos 9 + i cos 36» + 4 sin 9 Qy /l + cos26» ) 

4?r V 1 +cos 26» y - 3 cos 6»o +cos 36»o +4 sin e V l+cos 26»o 

f| n |- -(l+cos20o-sin20o+2Vl+cos20 o ) j _ j^j (1+cos 20Q-sin 20q+2V1+cos 20 o ) j j 

V"3 cos 2 o +cos O cos 30„+2 sin 20 o V^+cos W a/^3 COS 2 0()+COS (?o COS 30 o +2 sin 20 o Vl+cos 20 o 

+i sec 3 / 2 0o^/2 sin 2 20 o - (3 cos#o — cos 3^o ) 2 

[ln[ »(l+cos2g( ) -sin2g()-2Vl+cos26>o) j _ ^ -t(l+cos ggp-sin 20q-2V1+cos 20 (1 ) jj 

1 y/3 cos 2 gp-cos 0p cos 3e +2 si n 20 o V1+cqs~20q~ J 1 Vj} COS 2 00 — cos So cos 30o+2 sin 20 o %/I +cos 20o 

+2-\/2i VI + cos 26*0^-3 cos O + cos 36» + 4 sin 9 y/l + cos 26>i tan 2 6 }. 

(4-22) 



4.2 Asymptotic entanglement from nonlocal initial conditions 

Asymptotic entanglement levels for nonlocal initial conditions are reported for the first time in 
the context of the homogeneous Hadamard walk on the line. In order to show that the asymp- 
totic entanglement level is strongly dependent on whether the initial condition is localized or 
delocalized in position space. In this section we apply our inhomogeneous walk to consider in 
detail the case of initial conditions in the position subspace H p spanned by | ± 1), as 

|*(0)> = ^> 8 M. (4-23) 

Whit the projection on k-space this non localized state is simply \ip(k, 0)) = cos(k)(\L) +i\R)). 
The dependence on the initial conditions is contained in the coin operator. The time evolution 
of the motion in phase space for the initial states defined in the previous equation can be 
expressed by 

\M, 2t)) = (HqH^IM, 0)) for even times 

(4-24) 

\4>(k, 2t + 1)) = Hi(H Hiy\tp(k, 0)) for odd times. 
One can check that the eigenvalues of HqH\ are equivalent with Eq. (4-10), exactly and the 
normalized eigenvectors are of the form Eq. (4-11) with components 

u'(k) = sic e 2lfe - cis 
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Figure 3: (Colour online) Asymptotic of entanglement Se as a function of #0 and Q\ for 
ao(k) = ibo(k) = cos(fc) in even times. 

v(k) = —icqCi sin2A; 

w(k) = i\Jl — (c Ci cos 2k + s Si) 2 . 

As in previous subsection relate in details, the matrix elements of the reduced density operator 
after even times is obtained as 



^_ dku.ifu\\2f\F(k)\ 2 1 |G(fc)| 2 - 



+ 



A'-r 



= £*£{«'(*)«•(*)( 



l^(fc)l 2 j_ |G(fc)l S 



AT, 



— + ^^)+«'(A>*(£;)( 



Ar 



)}• 



(4-25) 



Subsequently, after odd times the relevant quantities for the entropy of entanglement are 

a = Z n §{(c?K(A;)| 2 + c lSl {u'{k)v*{k) + u'*(k)v(k)) + sl(\v(k)\ 2 + \ w (k)\ 2 ))(£gf + ^ 
+(c l8l (u'(k)w*(k) + u>*(k)w(k)) - 2s\v(k)w(k))(^f - ^f)}, 



P = J% ie 2tk {(cMW(k)\* - \v(k)f - \w(k)\ 2 ) - c\u\k)v*{k) + s \u'*{k)v{k)){^ + ^f-) 



+ (-c 2 u'(k)w*(k) + s\u'*{k)w{k) + 2c lSl v(k)w(k))( 1 -^ - ^gfi}. 

(4-26) 

From these two above equations, one readily obtains the independent elements of p c for 
this kind of initial conditions. Fig. 3 (A and B) shows the entropy of entanglement as a 
function of #0 and 9\ after even time steps and Fig. 4 (A and B) shows it after odd time 
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Figure 4: (Colour online) Asymptotic of entanglement Se as a function of 9q and 9± for 
ao(k) = ibo(k) = cos(fc) in odd times. 

steps. In comparison with local initial condition, the first aspect of these graphs to high- 
light is that, asymptotic entanglement in odd times versus both variations of coin opera- 
tor is symmetric and this distinguishing quality is free from type of initial conditions. The 
second aspect is that, S^ ter even Ume steps ^ vanishes which 9q tends to become |, also there 

ig Rafter even time steps) = QMm ^ ^ p ig 3) The origin of &U the aW factg 

/0 1 

can be clarified by using the analytical expressions. Let us begin by Hi = I and 

\1 

cos 9 sin 6*o \ 
H = ( , the wave functions of QW are 

sin 9 — cos 6> / 



|V>(2t)> 



£0! 



{(e- ite ° + (-l)V* e °)(|l) <g> \L) + i\ - 1) ® \R)) 



+( e -*«b _ _ 3) ® |L) + i|3) ® + - 1) ® |L) + z|l) ® |f2», 



[V(2* + 1)) = &-{(e- ite ° + - 2) <g> |L) + |2) ® \R)) 

+( e -«flb _ (_i)*e^o)( z | 2 ) g, |L) + | - 2) <g> |i?))} + if e -«*>|0) ® (i|L> + 



(4-27) 



Density operator after some manipulation, can be expressed as 



Pc (2t) 



and p c (2t + 1] 



1 i 

2 4 
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in the long-time limit, the contribution of the time dependent terms in these p c s vanishes and 
we have 

p c (after enen steps) = I J and p c (after odd steps) = I 

V I 2 / \ T 2 

From Eq. (3-6), the exact eigenvalues of the both operators are r± — | and r 2 = \, these 
eigenvalues yield the asymptotic value for the entropy of entanglement 

Se = -(I log 2 $ + ^ Iog 2 (^)) = 0.811278. 
If H replace by the shift coin operator, then the position states are 

|V(2t)> = t=fe ite i(\ - 1) + |1)) ® (|L) + 
|^(2* + 1)) = t^e^+D^lO) ® (\L)+i\R)). 

Walk is bounded in this cases, overtly, the reduce density operators are given by 

Pc (2t) =[*. M and p c (2t + 1) = I _ 2 . J 

V ~2~ 2 / V T 2 

we observe that entropy of entanglement is zero and 0.811278 in the even and odd times 
respectively. 



(4-28) 



5 Conclusion and discuss 

We have investigated an inhomogeneous walk with two-period, to be used in quantum random 
walk in one dimension latices. It determined by two orthogonal matrices and included two free 
parameters that together provided many conditions under which a measurement performed on 
the coin state yielded the value of entanglement on the resulting position quantum state. We 
have studied the problem analytically for all values of two free parameters of coin with diverse 
initial conditions for large number of QW steps. To summarise the results for this walk, we 
find the various behaviours of the entanglement are governed as follows. We demonstrated 
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that how the asymptotic value depend on coin parameters and initial condition. Moreover it 
was different for odd or even positions and when on the coin parameters take the precise value 
| walk altered to a bounded motion. Two striking characteristics of this walk was that, when 
motion began with local initial condition, walk regulated by the first coin operator that act 
on the state of particle in all even tims and S^ ter odd Ume steps ^ was symmetric ratio 6 and 9i 
for all type of initial conditions. 



Appendix 

The Riemann-Lebesgue lemma 

This little note is devoted to a proof of the Riemann-Lebesgue lemma. We use the following 
notation for the nth Fourier coefficient of a 27r-periodic function /: 

r7r dk ~, 



fin) = £ ^/(*)e** 



Lemma: Assume that f(k) is 27r-periodic, bounded and integrable. Then f(n) — > when 
n — > ±00. 

Proof: We shall prove this only for real valued functions. If f{k) is complex valued, the result 
will follow from the result applied to the real and imaginary parts of f{k) separately. First, 
we prove the result for an extremely special case: Namely, a single step, which is a function of 
the form 

'l a + 2k-K < k <b + 2kir, keZ 

otherwise 
where a < b and b — a < 2n. Then 



s(k) = < 



/ n f bdk mk e mb -e ma n 
sin) = / — e = > as n ^ ±00, 

since the numerator is bounded and the denominator goes to infinity. Second, since any step 
function is a linear combination of a finite number of single steps, the same result holds for 
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step functions. 

Finally, now assume that f(k) is integrable, and pick any e > 0. It follows - practically direct 
from the definition of integrability - that there exists a step function s(k) with 

From this we get 

_ w -(/(*) - ~s(k))e ink \ < J n ^\f(k) - m\ < e, 

as well. We have shown that s(n) — > 0, so there is some N so that |n| > implies s(n) < e. 
Whenever |n| > N, then 

< l/M-sMI + NMI <e + e = 2e, 

which finishes the proof. 

Notice that the Riemann-Lebesgue lemma says nothing about how fast f(n) goes to zero. 
With just a bit more of a regularity assumption on f(k), we can show that f(n) behaves 
roughly like 1/n or better. This is easy if f(k) is continuous and piecewise smooth, as is 
seen from the identity f'(n) = i nf(n), which arises from partial integration. Applying the 
RiemannLebesgue lemma to f'(k) we conclude that f(n) is 1/n times something that goes to 
zero, so f(n) — > faster than 1/n, [54]. 
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